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The thermal reactivity of HCN and NH 3 in interstellar ice analogues. 
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ABSTRACT 

HCN is a molecule central to interstellar chemistry, since it is the simplest molecule contain- 
ing a carbon-nitrogen bond and its solid state chemistry is rich. The aim of this work was to 
study the NH3 + HCN — > NH^CN thermal reaction in interstellar ice analogues. Laboratory 
experiments based on Fourier transform infrared spectroscopy and mass spectrometry were 
performed to characterise the NH^CN reaction product and its formation kinetics. This reac- 
tion is purely thermal and can occur at low temperatures in interstellar ices without requiring 
non-thermal processing by photons, electrons or cosmic rays. The reaction rate constant has 



a temperature dependence of k(T) - 0.016 



+0.010 

0.006 



exp( 



-2.7±0.4kJmor 



) when NH3 is much 



more abundant than HCN. When both reactants are diluted in water ice, the reaction is slowed 
down. We have estimated the CN ion band strength to be A CN - = 1 .8 + 1 .5 x 10~ 17 cm molec -1 
at both 20 K and 140 K. NHtCN exhibits zeroth-order multilayer desorption kinetics with a 



rate of k^ e& (T) = 10 28 molecules cm 1 s 1 exp( 



-38.0±1.4kJmor 



). The NH 3 + HCN 



NH+CN 

thermal reaction is of primary importance because ( i) it decreases the amount of HCN avail- 
able to be hydrogenated into CH2NH, ( ii) the NH^ and CN ions react with species such as 
H2CO, or CH2NH to form complex molecules, and ( Hi) NH4 CN is a reservoir of NH3 and 
HCN, which can be made available to a high temperature chemistry. 
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1 INTRODUCTION 

The HCN molecule plays a fundamental role in interstellar chem- 
istry as it is the simplest molecular species containing both C and 
N atoms (while the radical CN is the simplest interstellar species 
composed of C and N); reactions involving HCN can produce 
key species inter mediate to the forma tion of amino acids, such as 
aminoacetonitrile dDanger et al.ll201 lh . 

HCN has been observed in the gas phase in various en- 



vironments inclu ding molecular clouds dSnvder & Bun" 



IClarketalJ[l974h . circ umstellar enve lopes dBieging et 



Llhl 

al. 



1971; 

Il984l) . 



the comae of comets dlpetalj[l99(il) . the atmospheres of Nep- 
tune dRosenqvist et alj 19921) and Titan dTokunaga et al.ll98ll). and 
was r ecently discovered in the Murchison meteorite dPizzarelld 
l2012h . Typical abundance s with respect to H2 are 10 9 - 10 8 
towar ds molecular clouds dHirota et alJI 1998b and massive proto- 
stars dSchreveretalJll997h. rising to 10 6 in the hot envelope pre- 
cedin g a hot core dLahuis & van Dish oeck 2000; Bo onman et al.l 
l200ll) due to a combination of gas phase chemistry and the sub- 
limation of ices from dust grain surfaces. As such, HCN has been 
shown to be a linear tracer of ongoing star formation and is used 
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to probe dense gas in star forming regions dHirota et al.lll998h . The 
gas phase ratios between isomers such as HNC and HCN indicate 
the extent of photochemical isomerisation that has occurred in as- 
trophysical environments. The HNC/HCN ratio is highly depen- 
dent on te mperature, and has b e en observed to be ~ 1 in pre stel- 
lar cores dPadovani etai]|201ll: [Herbst et al.ll2OO0h . 0.6 - 0.9 in 
infrared dark clouds (Vasvunina et al 2011), and less than 0.3 in 
Neptune's atmosphere dMoreno et al. |201ll) . The gas phase ratio 
CN/HCN can change by orders of magnitude between the edge 
of a molecular cloud and the densest r egions of the core, as it i s 
sensitive to photodissociation channels dBoger & Sternberg 20051) . 
A greater knowledge of the grain surface chemistry network be- 
tween CN a nd HCN is needed to understand CN/HCN ratios in 
dense cores dChapillon et alj2012h . CN is the first diatomic anion 
to have been positively identified in an astrophysical environment, 
having been obser ved recently in the ga s phase in the C-star enve- 
lope IRC +10216 dAgundez et a"i]|2010l) . 

One environment where HCN is of particular interest is in 
comets and their comae. The CN moiety (par t) is key in the pro- 
duction of amino acids (Ferris & Hagan[[T984h , recent ly observed 
in com ets for the first time by the Stardust spacecraft (iE lsila et"al1 
2009). Both HCN and CN have been observed extensively in 
the gas phase on lines of sight towards cometary comae (e.g. 
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Huebneret aril 19741; iBockelee-Morvan et al]| 19871 ; IWoodnev et~ai1 
2002). HCN is believed to be the major source of the CN radi- 
cal dPaganini et al1l2010l) . although other sources, such as refrac- 
tory organic compounds, are necessary in order to produce the 
observed CN abundances jFravetalJ|2003) . In the solid phase, 
HCN is unlikely to exist solely in its monomeric form, but also 
in the form of polymers and complex species wh ich act as HCN 
reservoirs, such as hexa methyltetramine (HMT, iBernstein et al.l 
1991 IVinogradoff etalj|201ll) or poly(methylene-imine) (PMI, 
Danger et alj201 11) . Solid phase HCN has yet to be positively iden- 
tified in comets, but a tentative detection of solid pha se HCN on Tri- 
ton w as made recently using the AKARI telescope iBurgd orf et al.l 

HouJ). 

Experimentally, solid phase HCN has been produced under as- 
trophysically relevant conditi ons upon irradiation o f CH 4 :N 2 mix- 
tures with 5 keV electrons djamieson et al.l [2009t> . HCN can re- 
act as an electrophile, as in the reaction of HCN with H atoms i n 
the solid phase which produces methylamine dTheule et al.l2011ah . 
or as a nucleophile, as in its reaction with H2C O to form hy- 
droxyacetonitrile HOCHiCN dPanger et al.ll2012l) . The participa- 
tion of HCN in a low temperature solid phase Strecker synthe- 
sis (a series of reactions to synthesise an amino acid) yields 
aminoacetonitrile, NH2CH2CN, via t hermal processing of an ice 
mixture of CH 2 NH, NH 3 and HCN dPanger et alj|201 ll) . In addi- 
tion, the irradiation of HCN in ice mixtures and aqueous solutions 
with photons or charged pa rticles has been seen to yield comple x 
ions dGerakines et al.ll2004h . amin o acids dBernstein et al]|2002l) . 
amides, carboxylic acids and bases dColin-Garcfa et alj2009l) ! 

HCN is also a weak acid, therefore it could react with bases 
to produce salts. Previously, acid-base reactions have been shown 
to yield salts under astrophysically relevant conditions. The re- 
action HNCO + NIL produces NHIOCN^ dPemvk et alj 1 19981 ; 
Rauni er et al. I l2003al ; van Broekhuizen et al. 2004) with an activa- 
tion energy of 0.4 kj moL 1 dMispelaer et alfeoTlE ; HNCO + H 2 
produces H30 + OCN~ (Raunier et al. 2003b) with an activation en- 
ergy of 26 kj mor' (The ule et al .l2011al) . which then further reacts 
to form the isomeri sation product HOC N with an activation en- 
ergy of 36 kJ mor 1 dTheule et al.ll201 lbl) ; HCOOH + NH 3 yields 
NHJHCOCT dSchutteetalJll999h . Such species have relatively 
high sublimation temperatures (greater than ~ 150 K) compared 
with the simple molecules such as CO2, H2O and CO which com- 
pose the bulk of th e icy mantles on dust grains dRaunier et alj2004l ; 
iNoble et al.l2012ah . As such, it is likely that salts will remain on the 
grain surface following heating processes that sublime the majority 
of the ice mantle, ultimately formi ng part of the refra ctory organic 
residue observed on comet grains dGlavin et al.ll2008h . 

In this paper, we report the results of a study of the purely 
thermal reaction of the weak acid HCN with the base NH 3 : 



Table 1. Experiments performed in the current work. 



HCN + NH, NH+CN" 



(1) 



NH3 has bee n observed in ast r ophysical ices at abu ndances 
of 2 - 15 % H2O dLacv et alJI 19981 ; iBottinelli et ai]|2010l) . and is 
therefore a key basic, nucleophilic species available to react with 
acidic, electrophilic molecules in grain ice mantles. 

The 6.0 and 6.85 yum bands are observed in ices towards many 
astrophysical environments, including y oung stellar objects (YSOs, 
Keane et al. 2001; Booge rt et al. 2008) and quiescent molecular 
clouds 1 Boogert et al J201 ll) . but have yet to be fully characterised. 
NH^ is considered to b e a likely candidate to ac count for some of 
the feature at 6.85 yum dSchutte & Kha nna 2003). It has been sug- 
gested that most of the 6.0 jim excess, and at least some of the 



ID 



Molecules 



Ratio at deposition 









HCN 


NH3 


H 2 


C0 2 


CO 


i 


pure HCN 




1 














ii 


pure NH3 







1 











iii 


HCN & NH 3 




1 


20 











iv 


HCN & NH 3 




1.5 


1 











V 


HCN & H 2 




1 





18 








vi 


HCN, NH 3 & H 2 




1 


3 


30 








via 


HCN, NH 3 , H 2 & 


C0 2 


2.5 


6 


17 


1 





viii 


HCN, NH 3 , H 2 & 


CO 


2 


1 


2.5 





2 



' Ratios calculated from the experimental spectra. 



6.85 u rn band, is related to highly processed ices dGibb & Whittetl 
2002). However, given that the bands are observed in quiescent re- 
gions, some low-energy formation route to their carriers must be 
possible. Laboratory studies have shown that NH4 can be produced 
via acid-base reactio ns in ices at 10 K with no irradiation (e.g. 
iRaunieret al]|2003id) . 

Thus, reaction dTJ is important as it involves the abundant ice 
species NH 3 , produces NH^, and also involves HCN, the simplest 
molecule containing a CN bo nd, which could be a prerequisite for 
the formation of amino acids dDanger et alj|201 lh . 

In this work, the NH+CN~ product of the NH 3 + HCN 
reaction is spectroscopically characterised using Fourier 
transform infrared (FTIR) spectroscopy. The kinetics of the 
NH 3 + HCN reaction are investigated using isothermal ex- 
periments employing FTIR spectroscopy. The temperature 
dependence of the reaction rate is determined to be k(T) = 
0.016!°;,™ s-' exp( - 1J±0A R f mor ' ). The band strength of CN~ is 
determined to be A CN - = 1.8 ± 1.5 x 10~ 17 cmmolec -1 . Tem- 
perature programmed desorption (TPP) experiments using mass 
spectrometry allows calculation of the desorption rate as k<,\ cs (T) 



10 molecules cm 



exp( 



-38.0+1.4 kJ moP 



) for a zeroth-order 



multilayer desorption. NHJCN is therefore more refractory than 
H2O, and thus can participate in a warm, water-free chemistry. 



2 EXPERIMENTAL 

Experiments were perfo rmed using the RING experimental set-up, 
as described elsewhere dTheule et al.ll201 lbh . Briefly, RING con- 
sists of a gold-plated copper surface within a high vacuum chamber 
(a few 10~ 9 mbar). The temperature of the surface is controlled us- 
ing a Lakeshore Model 336 temperature controller, a closed-cycle 
helium cryostat (ARS Cryo, model PE-204 SB, 4 K cryogenera- 
tor), and a heating resistor; the temperature is measured using a 
PTGS 670 Silicon diode, with a 0. 1 K uncertainty. 

Molecular species in the form of room temperature gas mix- 
tures are dosed onto the gold surface by spraying at normal in- 
cidence via two injection lines. The infrared spectra of the con- 
densed mixtures are recorded by means of Fourier Transform Re- 
flection Absorption Infra Red Spectroscopy (FT-RAIRS) using a 
MCT detector in a Vertex 70 spectrometer. A typical spectrum has 
a 1 cirr 1 resolution and is averaged over a few tens of interfero- 
grams. In the gas phase, molecules in the high vacuum chamber are 
measured by means of mass spectroscopy using a Hiden HAL VII 
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Figure 1. The infrared absorption spectra of pure HCN and NH 3 deposited at 10 K: a) HCN at 10 K, b) HCN after heating to 100 K, c) NH 3 at 10 K, d) NH 3 
after heating to 100 K. 



RGA quadrupole mass spectrometer (QMS) with a 70 eV impact 
electronic ionisation source. 

In the experiments presented here, a series of gas-phase mix- 
tures of HCN, NH 3 , H 2 0, CO and C0 2 were prepared in the in- 
jection lines. Ammonia is commercially available as a 99.999 % 
pure gas from Air Liquide. The gas phase HCN monomer was syn- 
thesised via the thermal reaction of potassium cyanide, KCN, with 
an excess of stearic acid, CH 3 (CH 2 )[ 6 COOH, in a primary pumped 
vacuum line (Gerakines et al. 2003)- A small quantity of C0 2 was 
produced during HCN synthesis due to the decomposition of stearic 
acid. H 2 vapour was obtained from deionised water purified by 
several freeze-pump-thaw cycles carried out under vacuum. 

A complete list of experiments performed in this study is pre- 
sented in TableQ] In all experiments, NH 3 was dosed via one beam 
line, while the other species were dosed via a second line. Gas mix- 
tures were dosed onto the gold surface, held at 10 K. Except where 
noted in the text, the surface was subsequently subjected to a heat- 
ing ramp of 2 Kmin" 1 , and infrared spectra were recorded at inter- 
vals of 10 K. Molecules desorbing from the surface were monitored 
by mass spectroscopy. 



3 RESULTS 

3.1 Pure HCN and NH 3 

Figure [T] presents the infrared spectra of pure HCN (Experiment 
(i)) and pure NH 3 (Experiment (ii)) deposited at 10 K. The spectra 
have absorption bands at characterist ic frequencies, which are tabu- 
lated i n Table[2] The spectra of N H 3 (d'Hen decourt & Allamandolal 
1986t fSandford & Aliamandoialll993l) and HCN teernstein et al.1 
1997 ; lGerakines et al.1.2004 ) are well documented in the literature; 



NH 3 is characterised by four fundamental vibrational modes - v 4 
with frequencies 3212, 1064, 3375, and 1625 cm" 1 , respectively, 
while HCN displays vibrational modes vi - v 3 with frequencies 
3143, 822, and 2100 cirr 1 , respectively, as well as a prominent ab- 
sorption at 1617 cirr 1 assigned as the first overtone of the HCN 
bending mode. 

The pure HCN and NH 3 ices were heated from 10 - 150 K at 
a heating rate of 2 Kmin" 1 ; the spectra of NH 3 and HCN at 100 K 
are also presented in FigureQ] Between 70 - 80 K, crystallisation of 
NH 3 results in a rapid change in position and form of the four main 
bands: vi shifts to lower frequency (3206 cirr 1 ), v 2 and v 3 split 
into two peaks (1071 and 1093 cm" 1 , and 3378 and 3364 cm" 1 , 
respectively), and v 4 shifts to higher frequency (1649 cm" 1 ). These 
changes are visible in the spectrum of pure NH 3 at 100 K, shown in 
FigureQ] trace d. The crystallisation of pure HCN occurs gradually 
between 10 - 120 K, with no definitive change, unlike for NH 3 : v\ 
shifts to lower frequency (3139 cm" 1 ) at temperatures up to 60 K, 
before increasing to 3144 cm" 1 by 120 K; v 2 and v 3 both shift to 
higher frequency above 80 K (2099 an d 823 cm" 1 , respec tively, 
likely due to crystallisation of the HCN jDanger et~al]|201 ll) ); the 
overtone at 1617 cm" 1 shows similar temperature dependence as 
V\, shifting to 1620 cm" 1 by 60 K then returning to 1617 cm" 1 by 
120 K. Crystallisation is evinced by the difference between the low 
and high temperature spectra of FigureQ] 

The temperature programmed desorption of the pure species 
was measured by quadrupole mass spectrometry. The major frag- 
ments (relative intensity > 10 %) of NH 3 are the parent fragment 
(m/z 17, 100 %) and the loss of one hydrogen at m/z 16 (86 %); for 
HCN the major fragments are m/z 27 (100 %) and m/z 26 (17 %). 
The resulting QMS spectra of the TPD of pure NH 3 and pure HCN 
are plotted in the top panel of Figure|2] 

Thermal desorption is an activated process. The rate of des- 
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Table 2. Fundamental infrared band positions (cm '), for the species HCN, NH3, H2O, CO and CO2 in various mixtures (as defined in Table [T}. Unless 
otherwise stated, all values are for ices at 10 K. 



Vibration mode 


Assignment 


HCN 


NH 3 


HCN:NH 3 


HCN:H 2 


HCN:NH 3 :H 2 


HCN:NH 3 :H 2 0:X 








(i) 


(ii) 


(iii) 


(iv) 


(v) 


(vi) 


X = C0 2 (vii)° 


X = CO (viii) 


HCN/CN- 


v 3 


CN stretch 


2100 




2077 


2098 


2091 


2090 


2086 


2086 


NH 3 


v 3 


NH stretch 




3375 


3380 


3398 




3378 


3385 


3399 


NH 3 


V2 


umbrella 




1064 


1098 


1104 




1101 


1113 


1112 


C0 2 


V 3 


CO stretch 


2343* 


2342* 


2342* 


2343* 


2342* 


2342* 


2340 


2343* 


CO 


VI 


CO stretch 
















2138 


NHJ 


V4 


NH bend 






1479 


1471 






1482 


1498 


NHjClST r 


V4 


NH bend 






1435 


1435 




1437 


1437 


n/sf 


NHJCN- c 


v\ 


CN stretch 






2092 


2092 




2092 


2092 


n/a r 



" Deposition (and spectrum measured) at 15 K. * Present as a trace pollutant (produced during the HCN synthesis). Isolated on the surface at high 
temperature, after the desorption of NH 3 and HCN. The ice mixture containing CO (Experiment (viii)) was not heated to sufficiently high temperature to 

observe these bands. 



orption b y unit surfac e, r, can be expressed by the Wigner-Polanyi 
equation where the desorption rate constant fcdcs is de- 

scribed in terms of an Arrhenius law: 



dN 



:Ae~ E des /RT N" 



(2) 



where A is the pre-exponential factor, £des is the energy of des- 
orption of a molecule from the surface (Jmol -1 ), R is the gas con- 
stant (JKT 1 mol~'), T is the temperature of the surface (K), N is the 
number of adsorbed molecules on the surface (molecules cm -2 ), 
and n is the order of the reaction. The units of A depend on n: 
molecules '~" cm* 2 " s~ x . 

By rearranging Equation ® as: 



dN A -Ei irt .„ 

= - e des' TV" 

dT p 



(3) 



where jS=f =2K min 1 is the heating rate, the experimental data 
can be analysed. In order to calculate the energies of desorption of 
the isolated species HCN and NH 3 , zeroth order desorption kinetics 
are assumed, as is stan dard practise for the multilayer desorption of 
bulk material (see e.g. ICollings et al]|2004 iBurke & Browi]|20ld ; 
iNoble et alj2012ij) . In this case, A has units of molecules cm 2 s 1 . 

The energy of desorption, £ dcs was calculated using two meth- 
ods. Firstly (Method 1), assuming that the lattice vibrational fre- 
quency of the solid is 10 13 s _1 and the number of molecules 
in a monolayer is approximately 10 15 cnr 2 , the pre-exponential 
factor, A, can be fixed at a value of 10 28 molecules cnr 2 s -1 , 
and the desorption energy calculated as the only free parameter 
in the fit of Equat i on [3] t o the experimental data dCollings et al.l 
2003; Fuch s et al] 1 20061) . Using this method, the calculated 
desorption energies were: f^ds.HCN =30+1 kJmol -1 , and 
£ad S ,NH 3 =26+1 kJmor 1 . 

In the second me thod of fitting (Method 2), as proposed by 
lHasegawa etal] J 19921) the pre-exponential factor in Equation {3} 
was assumed to be a function of £ ads approximated by: 



A = jV ML .v = N ML 



2/v ML £ ads 



7T 2 M 



(4) 



where M is the mass of the adsorbate molecule, and 
Nml ~ 10 15 cm~ 2 . The advantage of this method is that the 



fit requires only one variable, £ des , rather tha n assuming or 
fitting the pre-exponential factor and fitting £des dAcharwa et al.l 
2007; INoble et alj l2012ah . Using the second method, the de- 
rived desorption energies were: ^ads.HCN = 28 ± 1 kJmol -1 , 
and £ a d s ,NH 3 = 25 ± 1 kJmol -1 . The pre-exponential factors 
were: A H cn =1.5 x 10 27 ± 5.2 x 10 23 molecules cnr 2 s -1 and 
A NH , = 1.7 x 10 27 ± 2.0 x 10 24 molecules cm" 2 s _1 . 

For both species, the values of the desorption energies cal- 
culated by the two m ethods fall within t he limits of uncer- 
tainty. As discussed in INoble et al.l d2012bl) , both methods are 
valid, since both pairs of parameters {A, E ies ] give the same 
value for the quantity of physical importance, i.e. the k Acs des- 
orption rate constant. We prefer to use the coupled solutions 
{A = 10 28 ± molecules cnr 2 s" 1 , E m 
and (A = 10 28 ± molecules cnr 2 s~', £ a 



30 



1 kJmol 



ads,NH 3 



26 ± 1 kJ moU 



This latter value is consi stent with the published valu e of 
25.6 ± 0.2 kJmoU 1 for NH 3 dSandford & Allamandolall 19931) . 



3.2 HCN and NH 3 mixtures 

The HCN and NH 3 gases were prepared in separate molecular 
beams and co-deposited on the gold surface held at 10 K, with an 
excess of NH3 in order to isolate the chemical reaction kinetics 
from the diffusion kinetics. Although this mixture (without H 2 0) is 
not directly astrophysically relevant, an excess of NH3 favours the 
complete reaction, gives HCN an homogeneous NH3 environment 
and thus produces reaction kinetics of (pseudo-) order one in HCN. 



3.2.1 Identification of the NHj + HCN reaction product 

The infrared spectra of Experiment (iii), a mixture of HCN and 
NH3, are presented in Figure [3] The bands attributed to HCN and 
NH3 are slightly shifted with respect to their values in pure HCN 
and NH 3 , as can be seen by comparing the values in Table [2] The 
presence of new bands in the spectrum at 10 K (c.f. Figure Q] and 
Table|2j also suggests that the species have undergone a reaction in 
the gas phase at 300 K during the deposition, before thermalisation 
to 10 K. It has previously been suggested that this reaction could 
also occur as a result o f the release of energy during condensation 
dGerakines et alj2004l) . 

The wide band at 1479 cnr 1 is attributed to the NHJ ion, 
which has been extensively studied in astrophysically relevant ices 
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Figure 2. The temperature programmed desorption spectra of a) pure NH3, 
b) pure HCN, c) NHJCN~. The experimental desorption spectra are plot- 
ted as solid lines (m/z 27 in black and m/z 17 in grey). Overplotted on these 
peaks are the best fitting zeroth order desorption kinetics for each molecular 
species; the fixed A method is plotted as a dotted line, while the parame- 
terised A method is plotted as a dashed line. The fitting approach is fully 
described in the text, where the two methods are denoted as Method 1 and 
Method 2, respectively. 



jRaunier et alJ l2003bL 120041 : iGerakines etai] |2004 iGalvez et al.1 
l20ld) . In addition to the broad band centred at 1479 cm 1 . which 
is assigned to the v 4 NH bending mode, NHJ exhibits a broad fea- 
ture in the 3000 cnr 1 region, which is the superposition of the v 3 
NH stretc hing mode, the 2v 4 overtone, and the y 2 + v 4 comb ina- 
tion mode dWagner & Horni Jl95Ct[ciutter & Thompsonll 19691) . In 
the 3000 cnr 1 region, at low temperatures (before NH 3 and HCN 
desorption), the NHJ contribution is difficult to deconvolve from 
the NH 3 contribution. At 150 K, when NH 3 and HCN have des- 
orbed from the surface, the IR spectrum of the product is isolated, 
and the NHJ feature is revealed as a three-peaked band with max- 
ima at 3168 (v 2 + v 4 ), 3014 (v 3 ), an d 2845 cm' 1 (2v 4 ). We ob- 
tain the same NH+CN~ IR s pectrum as lGerakines et al.l fc004l) and 
IClutter & Thompson! ( fl969l) . 

The band at 1479 cirr 1 grows and shifts redwards as the mix- 
ture of HCN and NH 3 is heated from 10 K to 150 K, as illustrated 
in the right hand panel of Figure [3] Upon the desorption of both 
reactants, HCN and NH 3 , between 1 10 and 150 K (as described in 
Figure[2} the band shifts to 1435 cirr 1 . This shift is due to both a 
change in environment and to crystallisation. 

The production of NHJ necessitates the presence of an anion 
in the ice. We positively identify the counter ion as CN~ from its 
IR band at 2092 cirr 1 (Figure [3] middle panel). As seen for the 
NH; band at 1479 cnr 1 , the HCN/CN~ band at 2077 cnT 1 shifts 
during heating (although it shifts bluewards), and between 1 10 and 
150 K the band narrows and shifts to 2092 cirr 1 due to a changing 
environment and crystallisation of the NH 4 CN~ product. 



We identify the product of the NH 3 + HCN reaction to be the 
NHJCN~ species with the following four arguments: 

• the identification of the NH^CN" absorption bands in the IR 
spectrum of the HCN:NH 3 mixture upon heating to 150 K (Fig- 
ure^ left hand panel), 

• in the TPD spectrum of the HCN:NH 3 mixture (shown for m/z 
27 in the lower panel of Figure[2] and discussed further in § I3.2.6K 
part of the product desorbs when dominant NH 3 desorbs around 
135 K. The mass spectrum of the pure product is obtained between 
140 K and 160 K, as shown in Figure[2] Mass features for the prod- 
uct are obtained at m/z 27 (100 %), 26 (17 %), 17 (20 %) and 16 
(15 %) at 150 K. The observed mass features of the product are a 
combination of the pure HCN and NH 3 desorption features. These 
features thus correspond to a proton transfer during the reverse re- 
action NH+ + CN~ — > NH 3 + HCN and the subsequent rapid des- 
orption of HCN and NH 3 at high temperature ( 140 - 160 K). Within 
the observational limitations of our experiments, no m/z 44 feature 
corresponding to the molecular ion NHJCN~ is observed, 

• the similarity of this system to the low tempe rature proton 
transfer reaction NH 3 + HNCO -> NH+ + OCN~ jDemvk et al.l 
1 19981 ; iRaunier et al]|2003al ; Ivan Broekhuizen et al.ll2004l) . 

• the previous literature evidence of the spontaneous reaction 
betw een HCN and NH 3 , represented in Equation l[T), upon deposi- 
tion dGerakines et al.ll2004klClutter & ThompsoiJ [l969 ) . 

We can therefore assign the product of the NH 3 + HCN reac- 
tion to the NH 4 CN~ species. 

A number of values can be derived from experiments on 
HCN:NH 3 mixtures: the band strength of the CN~ ion, A C n, the 
activation energy of the HCN + NH 3 reaction, and the desorption 
energy of the salt NH^CN" from the surface. 



3.2.2 The HCN bending mode at 848 cm' 1 . 

In a mixture of HCN:NH 3 , the only isolated HCN band (which 
does not overlap with bands of NH 3 or CN~) is the HCN bend at 
848 cirr 1 . Figure [4] (right panel) illustrates the effect of the envi- 
ronment on this band; spectra are plotted for Experiments (i) pure 
HCN, (iv) low NH 3 concentration, and (iii) high NH 3 concentration 
at 10 K. The band widens and shifts redwards with increasing NH 3 
concentration. The presence of the broad bending mode absorption 
feature at 800 cnr 1 in the spectrum of H 2 prevents the study of 
the HCN bending mode in mixtures containing H 2 0. In a sample 
of pure NHJCN~ at 20 K (Figure [4] trace d), prepared by mixing 
gaseous HCN and an excess of NH 3 in the injection lines before 
dosing, the HCN bending mode is not se en, as all HCN is co nverted 
to CN~ before dosing onto the surface jDanger et alj|201lt) . In all 
of the experiments presented here (without H 2 0) the HCN bending 
mode is clearly observed, indicating that in all experiments some 
HCN remains unreacted on the surface at 10 K. 

Conversely, the v 2 and v 3 bands of NH 3 both increase in fre- 
quency with the addition of higher concentrations of HCN. The 
v 2 band shifts from 1064 cnr 1 in pure NH 3 (Experiment (ii)) to 
1098 cnr 1 with a low HCN concentration (Experiment (iii)) to 
1104 cnr 1 with a higher HCN concentration (Experiment (iv)), 
while the v 2 band shifts from 3375 cnr 1 to 3380 cnr 1 to 3398 cnr 1 
in the same experiments. 

The absorption strength of the HCN bending mode in pure 
HCN at 10 K was calculated from the data in Experiment (i) to be 
j4«hcn = 1.9 ± 0.7 x 10~ 18 cmmolec -1 , using the area under the 
HCN bending mode and the HCN CN stretching mode peaks, and 
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Figure 3. IR spectra of Experiment (iii), a mixture of HCN and NH3 with an excess of NH3 . Left panel: Upon deposition at 10 K, and after heating to 150 K. 
Middle panel: evolution of the V3 CN stretching mode of HCN/CN~ as a function of temperature. Right panel: evolution of the V4 NH bending mode band of 
NH+ as a function of temperature. 
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Figure 4. A comparison of the CN stretching mode (left panel) and HCN 
bending mode (right panel) absorption bands at 10 K for a) pure HCN (Ex- 
periment (i)), b) HCN in a low concentration of NH3 (Experiment (iv)), and 
c) HCN in an excess of NH3 (Expe riment (iii)). Trace d is a spectrum of 
pure NHJCN" at 20 K taken from iDanger et alj|201 lh . The dotted lines 
mark the positions of the band maxima in Experiment (i), and are included 
to guide the eye. 



the ab sorption strength of the CN stretching mode ^Bernstein et al.l 
[l992l). 



3.2.3 The CN stretching mode at 2100 cm 

It is particularly interesting to note the effect of the environment on 
the CN stretch at 2100 cm" 1 , as the band is very complex, and the 
contributions of the species HCN and CN~ are largely irresolvable. 
The spectra of Experiments (i) pure HCN, (iv) low NH3 concentra- 
tion, and (iii) high NH 3 concentration at 10 K are plotted in Figure|4] 
(left panel), as for the HCN bending mode discussed in the previous 
section. The band at 2100 cm 1 behaves very similarly to the HCN 
bending mode at 848 cm 1 , widening and shifting redwards with in- 
creasing NH3 concentration. It is evident from Figure|4]that in the 
different experiments, different components contribute to the over- 
all profile of the band. By comparison to trace d in the same figure 
(pure NHJCN~), it can be inferred that an increasing concentration 
of NH3 results in an increased conversion of HCN to CN~ in our 
mixtures at 10 K. However, due to the complexities introduced by 
the effects on the band position and profile of concentration and, as 
previously discussed with respect to Figure |2j, temperature, it has 
not been possible to resolve the contributions of the species HCN 
and CN~ to the CN stretching mode in a reliable and reproducible 
manner across all of our spectra. 



3.2.4 Determination of the CN ion band strength. 

In order to determine the band strength of the CN~ ion, it is neces- 
sary to calculate the CN~ ion column density in our ice mixtures. 
Due to the difficulty of deconvolving the contributions of HCN and 
CN" at 2092 - 2099 ctrT 1 discussed above, it is not possible to 
evaluate directly the number of CN~ ions in the ice from its ab- 
sorption band at low temperature. At temperatures above 140 K, 
all HCN has desorbed from the surface, and therefore the fea- 
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ture at 2092 crrT 1 (4.78 jim) derives solely from absorption by 
the CN~ ion. The assumption is made that Equation [TJ holds and 
therefore the number of CIST ions is equivalent to the number of 
NHJ ions. Thus, using the area under the absorption feature at 
2092 cirT 1 and the number of NHJ ions in the ice (derived from 
ten HCN + NH 3 experiments) the average absorption strength was 
calculated as A C n , 140 k = 1.8 ± 1.5 x 10" 17 cm molec -1 . The uncer- 
tainty is calculated from the spread of values in the ten experiments, 
and is significant, reflecting the complexity of deriving such values 
from laboratory data. 

Using a spectrum of pu re NH^CN~ a t 20 K (shown 
in Figure [4] trace d) from iDanger et al.l ( l201ll) . we cal- 
culated the absorption strength of CN~ at 20 K as 
Acn-,20 k = 1.8 ± 0.1 x 10~ 17 cmmolec 1 . This result is identical to 
the value calculated at high temperature, suggesting that the CN~ 
band in pure NH|CN~ is not highly susceptible to temperature 
changes in the range of 20 - 140 K. 

The band strength of the CN~ ion is comparable to the band 
strength of most other interstellar molecules, unlike the OCN~ 
ion, which exhibits a high b and strength (1.3 x 10~ 16 cmmolec 1 
van Broekhuizen et al. 2004) that compensates for its low abun- 
dance and makes its detection possible in the IR spectra of inter- 
stellar ices. 

Further calculations were performed on data from Experi- 
ments (iii) and (iv) in an attempt to deconvolve the contributions of 
HCN and CN~ to the CN stretch absorption at 2100 cm" 1 . Firstly, 
the number of HCN molecules in Experiment (iv) was calculated 
from its absorption band at 848 cirT 1 (Figure[4] trace b), using the 
absorption strength calculated in § |3.2.2| (assuming that the magni- 
tude does not change with environment), to be ~ 9.3 x 10 17 . Sec- 
ondly, using the absorption strength calculated above for A C n-,20k 
at 2100 cm" 1 , it was possible to calculate the number of CN~ 
radicals in Experiment (iv), and thus infer the number of HCN 
molecules as ~ 2.5 x 10 18 , approximately 1.7 times higher than that 
calculated by the first method. The large uncertainty on this value 
most likely derives from the assumptions made in the calculations: 
that the absorption strengths of the CN~ and HCN CN stretching 
mode absorptions, and the HCN bending mode absorption do not 
change with enviro nment. As has been previously demonstrated by 
iBorget et alj d2012h . the CN stretch in nitriles is very sensitive to 
environment, with the absorption strength of CN in aminoacetoni- 
trile changing by a factor of two between the pure molecule and a 
mixture of aminoacetonitrile and H2O (1:3). Thus we conclude that, 
while we can accurately determine the absorption strength of CN~ 
at low to intermediate temperatures in a pure salt, it is not possible 
from these data to reliably deconvolve HCN and CN~ contributions 
to the CN stretching mode absorption at 2100 cm 1 . 

The implication of this for future observations is that in the 
infrared spectra of ices where both HCN and CN~ species are 
present, their bands are likely to be superimposed upon one an- 
other. Given the large difference in absorpti on strengths of the two 
species - 5.1 x 10~ 18 cmmolec -1 for HCN teernstein et al.ll 19971) 
and 1.8 x 10~ 17 cmmolec -1 for CN~ - it is important that the band 
not be treated as one species, but that efforts are made to deconvolve 
the individual contributions of HCN and CN~. Further experimen- 
tal studies are required to address this complex problem. 



3.2.5 Determination of the reaction rate constant 
determine 



In order to determine the activation energy of the 
NH 3 + HCN -» NHJCN" reacti on, isothermal kin etic exper- 
iments were performed, following iBossa et ail d2009l) . A series 
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Figure 5. Time evolution of the normalised abundances of the NH3 reactant 
and the NHJ product at different fixed temperatures. Top panel: NH3 evo- 
lution at 80 K (diamond) and 90 K (triangle). Bottom panel: NH+ evolution 
at 60 K (square) and 80 K (cross). Overplotted on the experimental data are 
the fits of EquationfJ] 



Table 3. Experimental rate constant, k, as a function of temperature. 



Temperature (K) 


k(x 10" 4 s" 1 ) 


60 


0.54 ±0.11 


70 


0.91 ± 0.27 


80 


2.64 ± 0.88 


90 


4.35 ± 1.00 


100 


4.11 ±0.93 


105 


7.09 ± 0.48 



of isothermal experiments was carried out in which mixtures of 
HCNNH3 with large excesses of NH3 (the ratios were typically 
1:15) were deposited at 15 K, heated as quickly as possible 
(~ 10 Kmirr 1 ) to a specific temperature, and held at that fixed 
temperature for a number of hours, in order to study the reaction 
kinetics. The temperatures chosen were 60, 70, 80, 90, 100, 
and 105 K. Above 105 K, desorption of the reactant molecules 
becomes important, and below 60 K the reaction proceeds so 
slowly that the experiments become prohibitively long to carry 
out. The abundances of NH 3 and NH^ in the ice mixtures were 
monitored via their characteristic absorption features in the IR 
spectra at 3375 and ~ 1440 cm" 1 , respectively. The time evolution 
of the species' molar fractions for four experiments are plotted in 
Figure [5] The other experimental data is omitted for clarity. It is 
evident from this figure that an increase in temperature increases 
the rate of the reaction, particularly in the first few hours of each 
experiment. 

As NH3 was in excess with respect to HCN, each HCN 
molecule was surrounded by a homogeneous NH3 environment. 
Thus, we can express the rates of destruction of NH 3 and forma- 
tion of NHJ with first-order reaction kinetics, in terms of the molar 
fraction of the reactants: 
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Figure 6. Experimentally determined ln(k) as a function of the inverse of 
the temperature, and the fit against an Arrhenius law (dot-dashed line). 



rf(HCN) 

dt 

,/(NH 3l 
rftNH^CN") 

dt 



-£(T).(HCN) 
-£(T).(HCN) 
ik(r).(HCN) 



(5) 



The first-order kinetics in HCN implies that the molar fraction 
of NH 3 is constant. Although not necessary in the first-order case, 
it is convenient to express the kinetic equations in terms of unitless 
molar fractions, which allows the calculation of a reaction rate con- 
stant k in s~ 1 . Resolving Equation[5]f or the isothermal experimental 
curves gives a reaction rate constant at a fixed temperature. The fits 
of four of the experiments are overplotted on the experimental data 
in Figure [5] The reaction rate constants, determined at the chosen 
temperatures by analysis of the NHJ band at ~ 1440 cirr 1 only, are 
displayed in Table [3] 

The temperature dependence of the reaction rate constant is 
plotted in Figure [6] in the form of ln(k) as a function of the in- 
verse of the temperature. Assuming k{T) follows an Arrhenius 
law, a least-squares straight line fit to the experimental data al- 
lows the calculation of the activation energy of the reaction HCN 
+ NH 3 . We derive an activation energy of 2.7 ± 0.4 kjmol 1 
with a pre-exponential factor of 0.016^ jj^ s" 1 . The uncertainty 
on the calculated values includes the dispersion of the results 
measured at 80 K and 90 K, where two experiments were per- 
formed at each temperature. The value for the activation en- 
ergy is consistent with v alues measured pr e viously for different 
thermal reaction syst ems jBossa et alj|2009l : I Theul e et al. 2011b; 
Mispe laer et al.ll2012h . Similarly, as observed previously, the pre- 
exponential factor is low, which may be due to the solvent cage 
effect in the solid solutions, i.e. the need for the molecules to 
orientate within the solid into favourable steric configuration be- 
fore reaction. Therefore the reaction rate constant has a k(T) = 
0.016+°;™ s" 1 exp( " 2 - 7±0 -^ m ° r ' ) temperature dependence within 
the small temperature interval (60 - 105 K) where we have been 
able to measure it. The question of the extrapolation to lower tem- 
peratures is critical to interstellar ice chemistry, but difficult to ad- 
dress experimentally. 



3.2.6 Determination of the NH^CN desorption rate constant 

A TPD spectrum of NHJCN" (m/z 27 and m/z 17) is shown in the 
lower panel of Figure|2] It is evident in this TPD trace that some of 
the product (or unreacted HCN and NH3) co-desorbs with the dom- 
inant species between ~ 120 - 140 K (similarly to the pure HCN 
in the upper panel) and that the NHJCN" desorbs at higher tem- 
perature. It is not possible from the experiments presented here to 
determine at what temperature NHJCN" first starts to desorb from 
the surface, as in none of our experiments was the NHJ CN" present 
on the surface with no HCN. Thus in all TPD traces, the signals of 
the desorption of HCN and NUJCN" overlap. The dominant mass 
fraction for NHJCN" is m/z 27, but as mentioned in § |3.2.1| above, 
other masses observed desorbing at this temperature include m/z 
26 (17 %), 17 (20 %), 16 (15 %), and 15 (2 %). Under current ex- 
perimental conditions there is no evidence of the desorption of m/z 
44 or 43, and thus it is concluded that NHJCN" desorbs from the 
surface entirely as HCN and NH 3 , following the reverse reaction 
NH+ + CN" -> NH 3 + HCN, rather than as NH+CN". 

The zeroth-order shape of the m/z 27 NHJCN" curve in Fig- 
ure [2] is less marked since only a few monolayers of NUJCN" are 
desorbing. Using the two methods described in § |3.1| above (fitting 
to the m/z 27 desorption peak only, and taking the average of val- 
ues calculated from five desorption spectra) the energy of desorp- 
tion of NUJ CN" was calculated as £ a ds,NH 4 CN = 38.0± 1.4 kj mor 1 



(Method 1) and E, 



35.4 ± 1.3 kJmor 1 (Method 2, with 



1.3 x 10 



,27 



ads,NH 4 CN 

6.0 x 10 23 molecules cirr 



). When consider- 



ing the energies derived for all three molecular species, the values 
calculated using the fixed A method are consistently 6 - 7 % higher 
than those calculated using the second method (A as a function of 

£ads)- 

3.3 HCN and NH 3 mixtures diluted in H 2 

In order to investigate the reaction of HCN and NH 3 under con- 
ditions relevant to ices in the interstellar medium, H 2 must be 
included in the mixtures. HiO is the most abundant solid phase in- 
terstellar molecule, with abundances of 10" 4 with respect to H 2 . 

3. 3. 1 HCN:H 2 O ice mixture 

Experiment (v) investigates the effect of dilution of HCN in a water 
ice at a ratio of 1:18. It is clear from the band positions tabulated 
in Table|2]that the peak position of the CN stretching absorption in 
HCN is shifted redward (with respect to that of pure HCN) when 
diluted in a H 2 ice (to 2091 cirr 1 ), although the shift is not as 
marked as that seen for the HCN/CN" band in a HCN:NH 3 mixture 
(2077 cm" 1 , Experiment (iii)). As expected, no reaction is seen for a 
mixture of HCN:H 2 0, neither at 10 K nor during heating to 180 K, 
where both HCN and H 2 have desorbed. 



3. 3. 2 HCN:NH 3 :H 2 O ice mixture 

Experiment (vi) is a mixture of HCN:NH 3 :H 2 with concentration 
ratio 1:3:30. This experiment investigates the HCN + NH 3 reaction 
when heavily diluted in H 2 0, as would be expected in interstellar 
ices. The IR spectrum of this mixture at 10 K is shown in Figure[7] 
trace a. 

The addition of H 2 shifts the HCN/CN" and NH 3 bands with 
respect to their positions in HCN:NH 3 mixtures. The HCN/CN" v 3 
band is observed at 2090 cm" 1 , while the NH 3 v 2 and v 3 bands are 
centred at 1101 and 3378 cm" 1 , respectively. 
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Figure 7. Left panel: Infrared spectra of a) a HCN:NH 3 :H 2 mixture (Experiment (vi), 10 K), b) a HCN:NH 3 :H 2 0:C0 2 mixture (Experiment (vii), 15 K), c) 
a HCN:NH3:H20:CO mixture (Experiment (viii), 10 K). Right panel: shift of the HCN/CN" band as a function of the ice mixture. The dotted line marks the 
position of the CN stretching mode absorption band in pure HCN. 



For this mixture, the NH^ band at 1460 cm is not observed 
as a distinct peak immediately after deposition at 10 K, unlike in 
mixtures of HCN:NH3 only. There is perhaps a trace amount of 
NHJCN" initially present in the mixture, but it is likely that the 
dilution of HCN and NH3 in a H 2 matrix prevents their direct 
reaction in the gas phase during the deposition. During the heat- 
ing of the ice mixture, a wide and shallow NHJ band, centered 
at ~ 1490 cm 1 develops in the spectra. At temperatures of 160 - 
170 K (when almost all of the H 2 has desorbed from the sur- 
face), this band shifts redwards to 1437 cm" 1 . The HCN/CN" band 
at 2090 cnT 1 shifts slightly redwards (6-7 cm" 1 , to around 2086 
cm" 1 ) during the heating of the ice mixture, due to the H 2 environ- 
ment, the conversion of HCN to CN", the crystallisation of CN", or 
a combination of these effects. Upon desorption of the H 2 0, when 
NHJCN" is the only species remaining on the surface, the peak of 
the band shifts to 2092 cm" 1 . 
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3. 3. 3 HCN:NH 3 :H 2 O: C0 2 and HCN:NH 3 :H 2 O: CO ice mixtures 

Further experiments were performed on astrophysically relevant 
ice mixtures (Experiments (vii) and (viii)) to investigate the effect 
of including C0 2 and CO in the ice on the position of the CN" 
absorption band. We wanted to investigate the possibility that the 
CN" ion could c ontribute to the band at 2041 c m" 1 (4.90 yum), as- 
signed to OCS dGibb et alJkOQ^lDartoislI 20051) , or to the so-called 
"XCN" band at 2165 cm" 1 (4.62 fjm) observed towards many as- 
trophysical environments, including high and low mass YSOs (e.g . 
Soifer et all 1 19791: i Gibb et al.l l200d : Ivan Broekhuizen et alj|2005l : 
Aikawa et al.ll2012h . It has been previously observed that the ice 



Figure 8. The normalised abundance of NHJ produced during the temper- 
ature ramps in Experiments (iii) and (vi). The reaction between HCN and 
NH3 without H 2 (Experiment (iii), diamonds) occurs more rapidly than 
the reaction in a H 2 matrix (Experiment (vi), triangles). Values are nor- 
malised to the abundance at 100 K, before the desorption of NH3 and HCN 
from the surface begins. 



absorption features of ionic species can be very sensitive to their 



chemical environment dSchutte & Khannall2003h and could shift by 
a few wavenumbers. 

The ice mixture containing C0 2 at 15 K is shown in Fig- 
ure [7] (left panel, trace b), while the ice mixture containing CO at 
10 K is shown in trace c of the same figure. The HCN/CN" and 
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NH 3 absorption bands are shifted by the presence of these two 
species in the ice, as tabulated in Table [2] Figure [7j (right panel) 
shows the spectral region containing the HCN/CN" stretching ab- 
sorption for all three astrophysically relevant ice mixtures studied. 
The traces are labelled as in Figure [7j (trace a) HCN:NH 3 :H 2 0, 
trace b) HCN:NH 3 :H 2 0:C0 2 , trace c) HCN:NH 3 :H 2 0:CO). The 
CN absorption feature of HCN/CN" is observed at 2100 cm -1 in 
pure HCN, at 2091 cm 1 when diluted into H 2 0, and at 2086 cirT 1 
in HCN:H 2 0:C0 2 and HCN:H 2 0:CO mixtures, is around 45 - 
60 cm" 1 from the 2041 cm" 1 band and around 65 - 80 cm" 1 from 
the 2165 cm" 1 band. 

By comparison of the absorption features in Figure|7] and the 
band positions in Table(2] it is clear that the addition of H 2 0, C0 2 
and CO to ice mixtures containing HCN and NH 3 is not sufficient 
to induce a significant shift in the position of the HCN or CN" fea- 
tures. The presence of CO and C0 2 in the mixtures shifts the peak 
to slightly lower frequency (2086 cm" 1 ). In the spectra presented 
here, it is not possible to deconvolve the contributions of HCN and 
CN" to the absorption feature centred at ~ 2090 cm" 1 . No absorp- 
tion features are observed in the region of 2165 - 2175 cm" 1 . 

We conclude that CN" can not contribute to either the 
2041 cm" 1 band, assigned to OCS, or to the "XCN" band. If present 
in ices, it would likely be observed at or very close to 2092 cm -1 . 

Regarding the kinetics of the reaction between HCN and NH 3 , 
it is evident that dilution in a H 2 matrix slows the formation of 
NH + CN", as illustrated in Figure [8] Quantification of the effect of 
diffusion on the kinetics of a thermal reaction is extremely compli- 
cated, and will be the subject of a forthcoming article (Mispelaer et 
al., in preparation). In this article, we have not attempted to further 
study the kinetics of the HCN + NH 3 reaction in mixtures contain- 
ing H 2 0. 



4 ASTROPHYSICAL IMPLICATIONS 

In low mass dark molecular clouds, the average abundance of gas 
phase HCN i s around 1-1.5 times higher than that of HNCO (see 
references in lVasvunina et al]|201ll) . and thus if OCN" is present 
in interstellar ices, it is reasonable to assume that HCN and CN" 
could equally be present. Gas phase HCN is present in molecu- 
lar clouds at abundances (with respect to H 2 ) of ~ 10" 7 - 10" 8 , 
while solid phase H 2 0, the most abundant solid phase species, is 
present at abundances of ~ 10" 4 . We can thus assume that the abun- 
dance of HCN in the solid phase is appro ximately 0.1 % H 2 Q . 
With values of 5.1 x 10" 18 cm molecule" 1 dBernstein et al.ll 19971) 
and 1.8 x 10" 17 cm molecule" 1 (this work) for the CN stretch of 
HCN and CN", respectively, the band strengths of these molecules 
are not high enough to compensate for their expected low inter- 
stellar abundances. Thus, given the sensitivity of current infrared 
telescopes, these two molecules are unlikely to be detected in ice 
mantles. An HCN detection has been reported in Titan, where at- 
mospheric c hemistry has resulted in a higher abundance of solid 
phase HCN dBurgdorf et alj20ld) . Moreover, neither the HCN nor 
the CN" band position in different ice environments matches the 
4.90 fim (2041 cm" 1 ) band currently assigned to OCS, nor can these 
species contribute to the 4.42 /jm "XCN" band. However, NH + CN" 
may contribute to the 6.85 fim band through its NH + bending ab- 
sorption band, as discussed above. 

The chemistry of HCN is of primary importance in interstel- 
lar ices, since HCN is the simplest molecular species with a CN 
bond and, as such, HCN is probably at the origin of some prebiotic 
molecules. An illustration of the chemical network centred upon 



HCN and NH + CN is shown in Figure|9] HCN can be produced by 
the hydrogenation of t he CN radical, and th en hydrogenated into 
CH 2 NH and CH 3 NH 2 dTheule et alj|2011al) . As the CN" moiety 
of the NH + CN" species is more nucleophilic than HCN, it could 
be the source of a rich ice chemistry. Its recent detection in the 
gas phase highlights its potential importance, and validates stud- 
ies in to the formation routes of interstellar ions dAgundez et al.l 
l20ld) . NH + CN" reacts with H 2 CO to form hydroxyacetonitrile 
HOCH 2 CN dDanger et al.l |2012|) . NH + CN" can also react with 
the HCN hydrogenation product CH 2 NH to give aminoacetonitrile 
NH 2 CH 2 CN, which can go on to form glycine according to the 
Strecker synthesis. However, the competing reaction of NH + CN" 
with CH 2 NH can al so form the poly(m ethylene-imine) CN-(CH 2 - 
NH)„-H (CN-PMI, lDangeretaUl201lh . The reaction product de- 
pends on the starting concentrations of NH + CN" and CH 2 NH; an 
excess of CH 2 NH promotes the polymerisation of CH 2 NH and the 
production of CN-PMI, while an excess of CN" promotes the pro- 
duction of aminoacetonitrile (and ultimately glycine). It is likely 
that CN" will react with other electrophiles in interstellar ices. The 
reaction HCN + NH 3 — > NH + CN" is the preliminary step neces- 
sary to the reactions discussed above and, as such, it is important to 
quantify the kinetics of this rate limiti ng reaction. Indeed, without 
H 2 0, HCN does not react with CH 2 dWoorJl200ll) . 

It is worth emphasising that the reactions in Figure|9]are ther- 
mally activated, meaning that they do not need any non-thermal 
process in order to occur at low temperature in astrophysical ices. 
The competition and the resulting branching ratios for these ther- 
mal reactions are entirely determined by the temperature depen- 
dence of their rate constants. As the reaction NH 3 + HCN is the 
first step in the network, its kinetics are limiting to the overall net- 
work. The measured barrier to this reaction (2.7 ± 0.4 kjmol" 1 ) 
is a little high c ompared to the form ation of the salt NH + OCN" 
(O^kJmol" 1 , [Mispelaer et"alll2012t) . but lower than the forma- 
tion of more comp lex species such as aminomethanol (4.5 kJ mol" 1 , 
iBossa etltHl 20091) . This simplified network is not complete; in or- 
der to fully explain the kinetics surrounding NH + CN", it is nec- 
essary to include all reactions, in particular those originating from 
NH 3 , such as the NH 3 + CQ 2 react ion, which can be either ther- 
mally act ivated teossa et al. I l2008h or activated by low-energy 
electrons dBertin et alj2009h . 

In addition to constraining its reaction network, it is impor- 
tant to understand how much HCN can be stored in the NH + CN" 
form, whether in molecular clouds, or in higher temperature en- 
vironments such as comets. When the HCN hydrogenation prod- 
uct CH 2 NH reacts with HCOOH, rather than with NH+CN", the 
complex molecules hexamethyltetramine (HMT, CeH^Nt) or a 
second form of poly(methylene-imine) (HCOO-(CH 2 -NH)„-H) are 
formed. In particular, HMT is known to decompose under VUV 
irradiation, producing HCN molecules and CN radicals, and there- 
fore is a possible extended source of the observed HC N and CN in 
comets dBernstein et all 19951 ; IVinogradoff et alj|201ll) . 

At higher temperature the reverse reaction 
NH+CN" NH 3 + HCN/HNC can occur, releasing NH 3 
and HCN (and possibly HNC) into the gas phase. This type of 
isomerisation has previ ously been shown f or the HNCO/HOCN 
system in water ice dTheule et alj 1201 lbl) . The formation of 
NH+CN" allows HCN to be stored within a less volatile species 
and made available to a high temperature chemistry, when the 
species fragments and eventually desorbs. 

As a follow up to studying its formation kinetics, it would be 
interesting to study the NH+CN" photodissociation under a VUV 
field. Does this sort of irradiation release the simple molecules 
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Figure 9. A schematic of the simplified chemical network surrounding 
the NH3 + HCN -> NHJCN~ thermal reaction. The final products are a) 
hydroxyacetonitrile (HOCH 2 CN), b) aminoacetonitrile (NH2CH2CN), c) 
poly(methylene-imine) (PMI, R-(CH 2 -NH)„-H, where R = HCOO orCN), 
and d) hexamethylenetetramine (HMT, C6H12N4). 



HCN and NH3, as upon heating the salt, or rather yield a more 
complex chemistry? Moreover, the HCN/HNC ratio produced by 
the reverse reaction could be measured in the gas phase using gas 
phase spectroscopy. This type of experiment could aid interpreta- 
tion of gas phase HCN/HNC ratios in interstellar environments. 

The inventory of the thermal reactions involving the most 
abundant solid state species must be continued, and the kinetics 
of the rate limiting reactions should be investigated, in order to 
achieve a complete network of interstellar grain reactions. The 
competition of this series of coupled reactions must be quantified 
as a function of temperature and of initial abundance ratios in order 
to understand interstellar ice chemistry. 



5 CONCLUSION 

In this work, we have shown experimentally that the purely ther- 
mal reaction between HCN and NH3 produces NHJCN" in inter- 
stellar ice analogues at low temperature. NHJCN~ has a charac- 
teristic CN~ absorption band at 2092 ctrr 1 (4.78 yum) that can be 
slightly redshifted by the ice environment. However, the environ- 
mental effect is not sufficient that the CN~ band can match the 
2041 ctrr 1 (4.90 pm) band observed in the IR spectra of interstel- 
lar ices. In addition, at low temperatures, the contributions of HCN 
and CN~ to the observed CN stretching mode absorption are diffi- 
cult to differentiate. We have estimated the CN~ ion band strength 
to be A CN - = 1-8 ± 1.5 x 10~ 17 cmmolec" 1 at 20 K and at 140 K; 
when coupled with its expected low abundance in interstellar ices, 
it is unlikely that solid phase CN~ could be detected by current IR 
telescopes. We have measured the kinetics of the pure HCN and 
NH 3 reaction with an excess of NH 3 , and we have derived the tem- 
perature dependence for the reaction rate constant to be k(T) = 
s -i exp( -2.7 ± o.4kJmor' y When the reactants are diluted 



nm £+0.010 „-l „.,„/ -2.7±0.4kJmol 
u ' ulu -0.006 



in water ice, the reaction is slowed down. We have also measured a 
desorption energy of £des,NH 4 CN = 38.0 ±1.4 kJmoL 1 , assuming a 
pre-exponential factor of 10 28 molecules cm -2 s . 
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